Astronomy & Astrophysics manuscript no. 9996 


© ESO 2008 


October 20, 2008 





MOST* discovers a multimode 6 Scuti star in a triple system: 

HD61199 

M. Hareter^**, O. Kochukhov^, H. Lehmann^ V. Tsymbal''"*, D. Huber\ P. Lenz\ W. W. Weiss\ J. M. Matthews^ S. 
Rucinski^, J. F. Rowe^, R. Kuschnig\ D. B. Guenther'', A. F. J. Moffat^ D. Sasselov'^, G. A. H. Walker^, and A. 

Scholtz^o 



Institut fiir Astronomie, Tiirkenschanzstrasse 17, 1 180 Vienna, Austria 

Department of Physics and Astronomy, Uppsala University, Box 515, 751 20 Uppsala, Sweden 
Thiiringer Landessternwarte Tautenburg, 07778 Tautenburg, Germany 
Tavrian National University, Dep. Astronomy, Simferopol, Ukraine 

Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, B.C., V6T IZl, Canada 
David Dunlap Observatory, Department of Astronomy, University of Toronto, Toronto, Ontario, L4C 4Y6, Canada 
Department of Astronomy and Physics, St. Mary's University, Halifax, Nova Scotia, NS B3H 3C3, Canada 
Departement de physique, Universite de Montreal, Montreal, Quebec, QC H3C 3J7, Canada 
Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts, MA 02138, USA 

University of Technology, Institute of Communications and Radio-Frequency Engineering, Gusshausstrasse 25/389, 1040 Vienna, 
Austria 



Preprint online version: October 20, 2008 



ABSTRACT 



Context. A field star, HD61199 (V « 8), simultaneously observed with Procyon by the MOST (Microvariability & Oscillations of 
STars) satellite in continuous runs of 34, 17, and 34 days in 2004, 2005, and 2007, was found to pulsate in 11 frequencies in the 
S Scuti range with amplitudes from 1.7 down to 0.09 mmag. The photometry also showed variations with a period of about four days. 
To investigate the nature of the longer period, 45 days of time-resolved spectroscopy was obtained at the Thiiringer Landessternwarte 
Tautenburg in 2004. The radial velocity measurements indicate that HD 61 199 is a triple system. 

Aims. A 6 Scuti pulsator with a rich eigenspectrum in a multiple system is promising for asteroseismology. Our objectives were to 
identify which of the stars in the system is the S Scuti variable and to obtain the orbital elements of the system and the fundamental 
parameters of the individual components, which are constrained by the pulsation frequencies of the S Scuti star. 
Methods. Classical Fourier techniques and least-squares multi-sinusoidal fits were applied to the MOST photometry to identify 
the pulsation frequencies. The groundbased spectroscopy was analysed with least-squares-deconvolution (LSD) techniques, and the 
orbital elements derived with the KOREL and ORBITX routines. Asteroseismic models were also generated. 

Results. The photometric and spectroscopic data are compatible with a triple system consisting of a close binary with an orbital 
period of 3.57 days and a 6 Scuti companion (HD 61 199 A) as the most luminous component. The d Scuti star is a rapid rotator with 
about vsin;= 130 kms"' and an upper mass limit of about 2.1 Mq. For the close binary components, we find they are of nearly equal 
mass, with lower mass limits of about 0.7 Mq. Comparisons to synthetic spectra indicate these stars have a late-F spectral type. The 
observed oscillation frequencies are compared to pulsation models to further constrain the evolutionary state and mass of HD 61 199 
A. The orbit frequency of the close binary corresponds to the difference of the two 6 Scuti frequencies with the highest amplitudes - 
a coincidence that is remarkable, but not explained. 

Key words. Stars: variables: S Scuti; Techniques: photometric, spectroscopic, radial velocities 



1. Introduction 

The 6 Scuti variables are stars of intermediate mass located in the 
classical instability strip near the main sequence, which pulsate 
in radial and nonradial acoustic (p-) modes driven by opacity ef- 
fects; see review by Breger (2000). Asteroseismic modelling of 
6 Scuti stars has been a challenge, even in cases where dozens 
of pulsation frequ encies have bee n identified - e.g ., FG Vir 
IBreger et al.l (I20ci5l) and HD 209775 lMatth"ewsl (l2007l) - because 

* Based on data from the MOST satellite, a Canadian Space Agency 
mission, jointly operated by Dynacon Inc., the University of Toronto 
Institute for Aerospace Studies and the University of British Columbia, 
with the assistance of the University of Vienna and on spectra taken 
with the Coude-echelle spectrograph attached to the 2-m telescope of 
the Thueringer Landessternwarte Tautenburg. 
** e-mail: bar eter@astro . uni vi e . ac . at 



the theoretical eigenspectra are so densely populated with fre- 
quencies. Constraining the fundamental parameters of a 5 Scuti 
star in a multiple system through its orbital elements is one pos- 
sible avenue to a unique asteroseismic model of the star 

HD61199 (V = 7.97) was observed by MOST, a Canadian 
space telescope dedicated to high-precision photometry of bright 
stars, in 2004, because it happened to lie in the field of MOST's 
first primary science target, Procyon (Matthews et al. 2004j). 
Little was known about HD61199 beyond i ts HD Catalogue 
classification as an A3 star iGatewood & Hani (|2006|) give a par- 
allax of 5.2+0.9 mas and likely luminosity class of V (see their 
Table 2, in which HD 61 199 corresponds to AO 1098). The 2004 
MOST photometry showed clear oscillations with periods in the 
6 Scuti range, plus an obvious variation with a periodicity near 4 
days. If the latter was due to rotational modulation, it would have 



2 M. Hareter et al.: MOST discovers a multimode 6 Scuti star in a triple system: HD 61199 

Table 1. Tetf and log g of HD 61 199 derived from published Table 2. Details of the MOST photometry of HD 61199. 

Stromgren colours. 

^Year 2004 2005 2007 





Moon 


Ribas et al. 


Balona 


Napiwotzki et al. 


7'eff 


7811 


7843 


7745 


7664 




3.94 


3.94 


3.95 


3.80 


spectral type 


A7 


A7 


A7 


A8 



given the rotation period of the star and helped for more clearly 
interpreting and modelling any rotational splitting in the pul- 
sation spectrum. Because of this possibility, groundbased spec- 
troscopy was obtained in 2004 from ThiiringerLandessternwarte 
Tautenburg to investigate the HD 61 199 system further. The star 
was again monitored by MOST in 2005 and 2007 as part of sub- 
sequent Procyon photometric runs. 

Photometric indices are available for HD61199, which 
can be employed to obtain an initial estimation of th e stel- 
lar pa r ameters . We h ave used TempLoggTNG (iKupka & Bruntil 
1200 It iKaised l2006h and published Strom gren indices b — 
y = 0.125, mi ^ 0.192 ,ci = 0.931 and //g = 2.816 
jHauck & Merminiodll998h for various calibrations to calculate 
values of T^g ranging from 7664 K to 78 1 1 K and log g from 
3.94 to 3.80 (see Table 1 ), based on different calibrations (iMoonI 
Il986t iRibas etanil997t lB"alonalll984t iNapiwotzki et aT]ll993h . 
These values were derived from the combined light of the triple 
system and imply a spectral type of A7 or A8, significantly later 
than the HD classification. The uncertainties in the individual 
temperature determinations are at least 200 K, and in the surface 
gravities, at least 0.3 dex. 

HD 61 199 is kinematically similar to DG Leo dFremat et al.l 

120051) . where two Am stars form a binary system and a third 
more distant component is a 5 Scuti star. This system is bright 
enough to decide which star is pulsating by means of ra- 
dial velocity variations, what would need larger telescopes for 
HD61199. 

The MOST photometry and data reduction is presented in 
Section|2] followed by the frequency analysis in Section|3] The 
spectroscopy and their analyses are described in Section|4] and 
the resulting fundamental parameters and orbital solution are 
presented. Asteroseismology is attempted in Section|5] where 
the observed pulsation frequencies are compared to models. 

2. MOST Photometry and Data Reduction 

The Canadian MOST (Microvariability & Oscillations of STars) 
microsatellite houses a 15-cm telescope feeding a CCD pho- 
tometer through a single broadband optical filter For de- 
scri ptions of the MOS T mission before an d after launch, 
see iWalker etH] ( l2003l) and iMatthewsl (|2004 . MOST's Sun- 
synchronous 820-km polar orbit (period = 101.4 min) allows 
uninterrupted observations of a target located in the satellite's 
continuous viewing zone for up to 8 weeks. MOST now can col- 
lect, simultaneously, three types of photometric data: (a) Fabry 
Imaging - in which the telescope entrance pupil is imaged as 
an annulus (diameter 40 pixels), for very bright stars ("Fabry 
Targets"), minimizing the influence of satellite pointing jitter 
and cosmic ray strikes on t he photometry, and enhancing the 
total signal (IMatthewsl |2004<) . (b) Dkect Imaging - in which a 
slightly defocused image of a star (a "Direct Imagin g Target") 
is projected directly onto the CCD dRowe et al.ll2006l) . (c) Guide 
Star photometry - in which photometry of up to about 30 stars 
used for MOST's startracking attitude control system is pro- 
cessed on board and downlinked to Earth. Since MOST's launch 



Length of dataset [d] 


33.57 


16.97 


34.25 


Integration time [s] 


0.9 


1 


24.8 


Sampling time [s] 


6 


5 


37 


# of raw data-points 


108 912 


49 978 


78 176 


# of data-points after reduction 


89 592 


47 079 


60249 


point-to-point scatter [mmag] 


4.9 


4.2 


2.6 



in June 2003 and its commissioning, the MOST team has im- 
proved the pointing performance of the satellite, improving the 
Direct Imaging photometric performance. The Guide Star pho- 
tometry mode was implemented during the mission as well, but 
was not available during the 2004 observations of HD 6 1 1 99 pre- 
sented here. 

HD 61 199 was observed simultaneously with Procyon from 
Jan. 8 to Feb. 9, 2004 and from Jan. 24 to Feb. 10, 2005 as a 
Direct Imaging Target, and from Jan. 4 to Feb. 12, 2007 as a 
Guide Star Target. The integration times were about 1 second in 
the first two runs, while the effective integration time was 24.8 
seconds in the 2007 run (in which short exposures were stacked 
on board). Table |2| lists details for the photometry. The pointing 
performance and other aspects of MOST photometry during the 
2007 Procyon were the best achieved to that point in the mission, 
resulting in significantly reduced point-to-point scatter. 

For the Direct Imaging photometry of H D61199 in 2004 
and 2005, a reduction routine developed by iHuber & ReegenI 
(2008) was u sed, simi l ar to t he decorrelation technique devel- 
oped by Reegen et alj HqqH) for MOST Fabry Imaging, but 
which also takes into account satellite pointing errors. Both tech- 
niques decorrelate the signal present in the target pixels with 
those in the background pixels. The 2007 Guide Star photom- 
etry of HD61199 is handled in a similar way, but for these 
data, there is no direct information about the background sky 
in each Guide Star subraster. The on-board Guide Star data pro- 
cessing is performed as follows: A Guide Star is sampled within 
a 20 X 20 pixel subraster. The subrasters are read out roughly 
every second and the average of the first and last columns of 
the subraster are used for an on-board background correction. 
This procedure is repeated until the final integration time of the 
Primary Science Target in the field is reached. Then the individ- 
ual background-corrected pixel intensities are integrated. Only 
this sum and the Guide Star exposure time are stored in the FITS 
header of the Primary Science Target data file which is down- 
loaded. Residual stray light due to scattered Earthshine remains 
in the Guide Star light curves sent to Earth. Fortunately, the mod- 
ulation of the stray light with the MOST orbital period is well 
known, and 15 Guide Stars were observed simultaneously with 
HD61199 allowing for a modified decorrelation technique us- 
ing non-variable stars in the field to correct for residual instru- 
mental and background variations. 

Figure[T]shows the reduced MOST Hght curves of HD 61 199 
for the three observing runs and their respective spectral win- 
dows are shown in Fig.|2| The low-amplitude alias sidelobes 
visible in the spectral window of the 2007 data, spaced by the 
MOST sateUite orbital frequency (164;uHza! 14.192t/-'), are 
due to our more conservative rejection of outliers in the Guide 
Star reduction. This procedure is producing more gaps when 
stray light is highest at certain phases in the MOST orbit than 
in the 2004 and 2005 Direct Imaging reductions. 

The decorrelation method reduced the amplitudes of 
orbitally-modulated background variations in the Guide Star 
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Fig. 1. MOST light curves of HD 61199 from 2004 (top panel), 
2005, and 2007 (bottom panel). The solid line shows a running 
average of 600 data points. 



.-50 -40 -30 -20 -10 10 20 30 40 50 



, .-50 :40 :30 :20 -10 10 20 30 40 50 



-50 -40 



-30 -20 -10 10 20 30 40 50 

Frequency [c/d] 



Fig. 2. Spectral windows for 2004, 2005 and 2007 from top to 
bottom. Aliasing is practically not present in the Direct Imaging 
data of 2004 and 2005. The more restrictive outlier rejection pol- 
icy in the Guide Star photometry of 2007 affects the spectral 
window more than the Direct Imaging reduction. 



data to about 8% of the raw values. For the final light curve of 
the 2007 run, a second-order trend (possibly due to long-term 
stellar variability or subtle instrumental effects) was subtracted. 
The subtraction of this polynomial fit affects neither the iden- 
tified frequencies nor their amplitudes, which were checked by 
comparing the discrete Fourier transform spectra before and af- 
ter subtraction. 



3. Frequency Analysis 

SigSpec dReegenl 12007 ) and Period04 (Lenz & B regeJ l2005h 
were adopted for the frequency analysis. SigSpec performs itera- 
tive prewhitening of the data with the most significant frequency 
found in each iteration. The iterations end if a specified thresh- 
old is reached ; the default value is 5.5 in spectral significance 
( lReegenll2007! ). corresponding to a signal-to-noise (S/N) ratio in 
amplitude of about 4. A peak identified with a significance of 
5.5 means that a similar peak has a probability of only 10"^^ of 
being produced by white noise (i.e., once in 3.2 x 10^ cases). 
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Fig. 3. Fourier amplitude spectra of the 2004 (top), 2005 (mid- 
dle) and 2007 (bottom) light curves. Power due to the stray light 
modulation with the MOST orbit is present, but with signifi- 
cantly reduced amplitudes compared to the raw data. The ver- 
tical dashed lines indicate the MOST orbital frequency, its har- 
monics and its 1 d ' sidelobes. (The MOST Sun-synchronous 
orbit carries it above nearly the same point on Earth after 1 day, 
resulting in a secondary modulation of the scattered Earthshine.) 



Fig. 13 shows the Fourier amplitude spectra of the individ- 
ual observing runs. Frequencies were searched in the range 
- 80t/"', down to a significance limit of 5.5. The final set 
of 13 adopted frequencies listed in Tab.[3] includes only those 
which could be detected in all three runs, where 11 of the 
frequencies are in the S Scuti domain, and the other two are 
close to 0.25 1/"'(/3 and fi). The uncertainties of the frequencies 
and amphtudes are c alculated according to the technique of 
iKallinger et all (l2008l) . Based on extensive simulations, those 
authors find that the frequency errors are overestimated when 
using the classical Fourier T^^ criterion. They propose instead: 



and cr(fl) - 



where cr(f) and cr(a) are the errors in frequency and am- 
plitude, respectively; T is the time base of the observations 
in days; and sig{a) denotes the significance of an individual 
frequency with amplitude a. 

The amplitudes a in Table[3]were calculated from the com- 
bined light of the system, which we know to be triple (see 
Sec.|4|. Therefore, the intrinsic amplitudes are actually higher 
by the inverse fraction of light of the pulsating component. We 
show in Sec.|4]that about 80% of the total flux can be attributed 
to the 6 Scuti star in the system. Hence, a scaling factor of 1 .25 is 
needed to obtain the intrinsic amplitude. In addition, the pulsa- 
tion amplitudes do show evidence of variability across the three 
runs, particularly for /4, as can be seen in Fig.|4] At this stage 
of the investigation, we can only speculate if this variability 
is due to unresolved beating frequencies or due to another ef- 
fect. In a recent study, Breger & Lenz investigated possible ex- 
planations for amp litude variability in the 6 Scuti star 44 Tau 
( iBreger & Len3l2008i) . 

The two close frequencies fg and /n in Fig.|4] are spaced 
0.063 d ' (for 2004 data). The frequency errors given for both 
frequencies in Table[3] are 0.005 and 0.009 respectively. Even 
the more conservative estimation, the T"' criterion gives a fre- 
quency resolution of 0.029 d '. Hence the frequencies are re- 
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Table 3. List of identified frequencies and amplitudes in the three MOST data sets. 



label 


frequency 


amplitude 


signif. 


frequency 


amplitude 


signif. 


frequency 


amplitude 


signif. 




2004 [rf-'] 


2004 [mmag] 


2004 


2005 [d-'] 


2005 [mmag] 


2005 


2007 [d-'] 


2007 [mmag] 


2007 


/i 


25.257 


±0.001 


1.72 ±0.05 


872 


25.254 


±0.003 


1.53 ±0.07 


534 


25.256 


±0.001 


1.46 ±0.05 


930 


/2 


24.974 


±0.002 


0.98 ±0.06 


308 


24.972 


±0.003 


1.19 ±0.06 


344 


24.973 


±0.001 


1.14 ±0.05 


614 


/3 


0.254 


±0.007 


0.84 ±0.06 


231 


0.241 


±0.015 


1.06 ±0.06 


287 


0.255 


±0.007 


0.98 ±0.04 


512 


/4 


22.625 


±0.002 


0.79 ±0.06 


206 


22.624 


±0.005 


0.78 ±0.06 


156 


22.624 


±0.001 


1.14 ±0.05 


635 


/5 


27.628 


±0.003 


0.62 ±0.05 


131 


27.627 


±0.005 


0.70 ±0.06 


130 


27.625 


±0.002 


0.70 ±0.04 


307 


/6 


23.582 


±0.003 


0.53 ±0.05 


99 


23.579 


±0.005 


0.67 ±0.06 


119 


23.582 


±0.002 


0.73 ±0.04 


327 


/7 


0.286 


±0.007 


0.42 ±0.05 


60 


0.275 


±0.017 


0.20 ±0.06 


12 


0.279 


±0.007 


0.65 ±0.04 


262 


/s 


18.697 


±0.007 


0.23 ±0.05 


19 


18.691 


±0.014 


0.16 ±0.04 


17 


18.691 


±0.004 


0.25 ±0.04 


44 


/9 


23.435 


±0.005 


0.35 ±0.05 


43 


23.441 


±0.015 


0.24 ±0.06 


15 


23.437 


±0.006 


0.18 ±0.04 


24 


/lO 


29.697 


±0.006 


0.27 ±0.05 


26 


29.689 


±0.016 


0.22 ±0.06 


13 


29.701 


±0.005 


0.24 ±0.04 


40 


/u 


22.329 


±0.010 


0.15 ±0.05 


8.6 


22.329 


±0.023 


0.16 ±0.06 


6.5 


22.326 


±0.007 


0.16 ±0.04 


18 


/l2 


23.372 


±0.009 


0.18 ±0.05 


12 


23.359 


±0.025 


0.14 ±0.06 


5.4 


23.373 


±0.012 


0.09 ±0.04 


6.4 


/l3 


24.128 


±0.009 


0.17 ±0.05 


11 


24.158 


±0.026 


0.13 ±0.06 


5.1 


24.130 


±0.009 


0.12 ±0.04 


10 
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Fig. 4. Frequencies in the 6 Scuti range detected in HD 61199. 
Uncertainties in amplitude (see text) are approximately the sizes 
of the symbols. 



solved. /i2 is close to an expected sidelobe with order 5 (i.e. 
2 X 14.192 - 5 c/"') of the first harmonic (28.384 d"') of the 
MOST orbit frequency. Sidelobes tend to decrease in amplitude 
with growing order. No sidelobe has been detected below the 
orbit harmonic of 28.384 d except for order 1 at 27.387 d 
Thus, /i2 has been adopted in the frequency list (Table[3]). 

The original 2005 light curve of HD61199 shows longer- 
term variability with a period of about 4 days which can be 
seen by eye (and is not present in the other stars observed in 
the same field). This timescale was investigated through Phase 
Dispersion Minimization (PDM, Stellingwerf 1978). The PDM 
diagrams of the 2004, 2005 and 2007 photometry are plotted 
in Fig.|5] For this analysis, all significant power at frequencies 
above 1 c/ ' has been prewhitened and the data have been av- 
eraged in 1-min bins. The most prominent feature in all three 
PDM diagrams occurs near a period of about 4 days, or more 
precisely, of 3.94 + 0.12 days in the 2004 and 2007 data and of 
4. 15 + 0.24 days in 2005. For those eiTo rs, the error estimation 
of 0.25r ' is adopted, which is given bv lKallinger et alj (|2008|) 
for the case of two closely spaced frequencies. The change in 
the period is within the errors. Prewhitening of the dominant pe- 
riod in each data set yields the PDM results given in the upper 
curves in each panel. This reveals clearly in the 2007 data, the 



1 

0.98 
0.96 



' '. ' 1. ' 1 ' ' 1 ' ' ' ' 1 




. 1 . 1 1 1 1 


11 1 11 11 _ 






















, 1 , < < < 1 






0.98 
0.96 



4 



j_IJ 

3.5 4 4.5 

Period [d] 



Fig. 5. Phase Dispersion Minimization (PDM) diagrams for the 
2004 (top), 2005 (middle) and 2007 (bottom) photometry. The 
lower curve in each panel is the PDM after prewhitening all sig- 
nificant frequencies higher than It/'. Each upper curve shows 
the PDM of the residuals after prewhitening by the dominant pe- 
riod in the lower curve. The thin vertical line marks the orbital 
period of the stellar system found spectroscopically (see Sec.|4]i. 



stellar system orbital period determined by spectroscopy (3.57 d, 
shown by the vertical line, see Sec.|4]i. The period may be present 
marginally in the PDM analysis of 2004 data, but is absent in the 
2005 PDM residuals. The SigSpec and Period04 routines find 
this latter period (fi in Table|3ll in all three data sets, within the 
expected frequency uncertainty. 



4. Spectroscopy 

Prompted by the longer period seen in the 2004 MOST pho- 
tometry, a spectroscopic follow-up campaign was carried out. A 
total of 29 high-resolution spectra (R-63 000, wavelength range 
from 4760 - 7360 A) spanning 45 days were obtained with the 
Coude-echelle spectrograph attached to the 2-m telescope of the 
Thiiringer Landessternwarte Tautenburg. The Signal to Noise 
(S/N) of the spectra is between 60 and 180 with a mean value 
of 130. The S/N was measured in the continuum around 5200 A. 

We analysed the dynamics of the HD61199 system by 
applying the leas t-squares deconvolution technique (LSD, 
iDonati et al.l 1 19971) to the spectra. Starting with preliminary 
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Table 4. Radial velocities for each component of the HD 61199 
system derived from the Tautenburg spectra. 



No. 


HJD 


Phase 






Va 




-2453000 




(kms-') 


(kms-i) 


(kms-') 


01 


69.2601 


0.048 


109.3±0.3 


-33.3±0.4 


46.8+1.1 


02 


69.2818 


0.054 


108.7±0.3 


-33.0+0.4 


47.3+1.1 


03 


69.3038 


0.060 


107.3±0.4 


-31.3±0.4 


48.1 + 1.4 


04 


70.2833 


0.334 


1.0±0.3 


77.6±0.4 


45.7+0.9 


05 


70.3049 


0.340 


-1.4±0.3 


80.3+0.4 


46.6+1.0 


06 


76.3648 


0.036 


110.8±0.7 


-36.6+1.0 


44.7+2.5 


07 


79.3111 


0.860 


85.1 ±0.4 


-10.3+0.5 


45.9+1.4 


08 


80.3180 


0.142 


85.7±0.5 


-11.4+0.6 


45.4+1.5 


09 


81.2756 


0.409 


-24.1±0.3 


103.6+0.4 


45.4+1.2 


10 


82.2736 


0.689 


9.8±0.3 


68.0+0.4 


43.7+0.9 


11 


86.3425 


0.827 


73.3±0.4 


1.8+0.5 


46.6+1.2 


12 


87.3380 


0.105 


96.9±0.4 


-23.3+0.5 


46.1+1.3 


13 


90.2719 


0.926 


105.1±0.6 


-28.8+0.4 


45.7+2.3 


14 


91.3869 


0.238 


43.3±1.0 


33.7+1.3 


45.5+1.2 


15 


91.4086 


0.244 


42.7±1.3 


35.1+1.1 


46.7+1.3 


16 


93.3190 


0.779 


51.7±0.4 


24.5+0.5 


45.1+1.0 


17 


93.3409 


0.785 


54.4±0.3 


22.1+0.4 


44.6+0.9 


18 


93.4236 


0.808 


64.3±0.3 


11.0+0.4 


45.0+1.0 


19 


95.3082 


0.335 


-0.1±0.3 


78.1+0.4 


46.5+0.9 


20 


96.3153 


0.617 


-17.3±0.3 


96.0+0.4 


47.0+1.1 


21 


97.3154 


0.897 


97.6±0.3 


-23.8+0.4 


46.6+1.1 


22 


98.3171 


0.177 


71.4±0.4 


2.6+0.5 


44.2+1.0 


23 


99.3164 


0.457 


-33.1±0.3 


113.0+0.4 


46.7+1.1 


24 


107.3302 


0.699 


14.6±0.5 


63.5+0.6 


45.2+1.5 


25 


110.2993 


0.529 


-34.5±0.3 


114.2+0.4 


46.7+1.1 


26 


111.2962 


0.808 


64.3±0.3 


10.7+0.4 


44.0+0.9 


27 


112.3031 


0.090 


100.5±0.4 


-27.7+0.4 


46.3+1.1 


28 


113.2990 


0.369 


-16.3±0.4 


91.6+0.5 


45.7+1.1 


29 


114.3325 


0.658 


-3.3± 0.3 


80.7+0.4 


45.4+1.0 



estimates for the atmospheric parameters, we extracted from 
the Vienna Atomic Line Database (VALD. iKupkaet aLlll999l: 
iRyabchikova et al.ll997LlPiskunov et al.ll995h 700 spectral lines 
with a residual intensity higher than 0.2 and used this list for the 
LSD line mask. Based on Table[Tl following model parameters 
have been adopted: T^ff - 7800 K, log g - 4.0, and solar chemical 
composition. The relative shape of the LSD profile is rather in- 
sensitive to line masks with different parameters. Thus, because 
we are using only relative quantities, no additional error is intro- 
duced by the difference of Teff of the components. Fig. |6] shows 
the LSD profiles of all 29 spectra with a typical S/N of about 
500. A broad - but shallow - line of a rapidly rotating star (to 
which we refer as component A) is clearly present, superposed 
on the narrow lines of an obvious binary system (referred to as 
components Ba and Bb). At some phases (e.g., spectrum #26) a 
fourth faint narrow component appears at the constant velocity 
of -26kms-'. At the phase of 0.926 (spectrum #13) it blends 
the lines of one binary component. For the spectra in which the 
narrow line component is present, the heliocentric velocity of 
the observer proj ected in the direction of HD 61199 varied in the 
range of -25 to -28 kms^^ Hence, this faint component may well 
be due to stray moonlight. The combined light of the system is 
dominated by the rapidly rotating star, which causes the broad- 
line component in the LSD profiles. This feature did not change 
significantly its position and its equivalent line width is much 
larger than those of the narrow line components. 

We fitted the LSD spectra with a simple model consisting 
of three superimposed Gaussian profiles, each broadened with 
a FWHM - 5 kms" ' . The following free model parameters were 
adjusted using a least-squares procedure: vsin/, radial velocity 



and intensity. The LSD profiles were fitted in a window between 
-140 and -1-235 kms"' and the resulting best-fit solution is shown 
in Fig.|6] The radial velocities measured from the spectral fits 
are listed in TableU] Errors of the LSD profiles were estimated 
from the discrepancy between the LSD model and observations. 
Uncertainties of the parameters determined from the fit to LSD 
profiles are formal errors calculated by a least-squares fitting 
routine. 



4.1. Determination of tine dose Binary Orbit 

To determine the orbital elements for the close binary evident 
in the spectra, two different approaches were used: (i) analysis 
of the radial velocities determined from the fits to LSD profiles 
using ORBITX (Tokovinin 1992) and (ii) disentangling the spec- 
trum with KOREL (Hadrava.l995h . 

Due to similar line intensities of the components Ba and 
Bb it is not straightforward to identify their contribution in 
the composite line profile. Hence a time-dependent observable, 
IVfia - Vfibl, was formed which varies with twice the orbital fre- 
quency of the close binary. Fig. [T] (top panel) shows the result- 
ing Fourier amplitude spectrum. The two highest peaks occur 
at frequencies of 0.4375 and 0.5616ii-', corresponding to peri- 
ods of 4.545 and 3.561 days, respectively. Phase plots with both 
possible periods are shown in the middle and bottom panels of 
Fig.|2] A period of 4.55 days leads to discontinuities in the phase 
diagram, while a period of 3.56 days produces a coherent RV 
curve. The correct identification of components Ba and Bb ap- 
pears to fail only for spectrum #13 where the the constant feature 
at -26kms-' interferes, and possibly in spectrum #15 where the 
LSD profiles of the narrow line components are blended. 

Second, KOREL was applied with the following setup: three 
stars forming two hierarchical orbits, with components Ba and 
Bb close together and A in a wide orbit. Spectra of the rapidly 
rotating B3V star HD 120315 were used to model the spectrum 
of telluric lines which were subtracted from the broad-lined 
component to obtain a pure stellar spectrum. Orbital elements 
were made free parameters in the fit and line strengths were 
allowed to vary with time. Since the line positions of compo- 
nent A did not change measurably during the period of observa- 
tions, we assumed the period of the wide orbit to lie in the range 
1 000 - 20 000 d. For 12 different long trial periods, we detected 
only marginal changes in the disentangled spectra and in the re- 
sulting orbital elements. 20 000 d is the upper limit to observe 
any changes in the KOREL solution. 

KOREL determines simultaneously orbital elements, line 
strengths of each component and disentangled spectra. Radial 
velocities are not determined from the input spectra but are cal- 
culated from the orbital solution. The KOREL code finds several 
local minima that all give zero RMS of the obtained RVs (shifts 
applied by KOREL to the individual spectra), i.e. the RVs are 
identical to those calculated from the simultaneously obtained 
orbital elements. We searched for the best solution nearest to that 
obtained with ORBITX, using the wavelength interval 4895 - 
5680 A. For an error estimation we then divided this interval into 
10 parts of equal lengths and applied KOREL to each of them. 
In this way we obtained 10 orbital solutions, mean and standard 
deviation, which we weighted by the total equivalent width in- 
cluded in each subspectrum. There was no significant difference 
between the values obtained from the full wavelength interval 
and the weighted means. 

The disentangled spectra and the combined spectrum are 
shown in Fig. [8] KOREL provides no information about the in- 
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Fig. 6. LSD profiles phased with the fi- 
nal derived orbital period of 3.5744(3) 
days. The red solid curve shows the fit 
of the triple- spectral model. The spec- 
trum numbers (in chronological order) 
are given on the left of each spec- 
trum, while orbital phase is given on 
the right. 



dividual continua of the stars, because the disentangled spectra 
are normalised to the common continuum of the observed com- 
posite spectrum. To renormalise the KOREL output spectra, one 
must either use information about the luminosity ratios of the 
stars or, as we did, scale the stellar spectra by comparison to 
synthetic spectra. KOREL has difficulties disentangling the H„ 
region, where the RV shifts are much smaller than the intrinsic 
line width. For this reason the disentangled spectra cannot be 
used for Teff estimation from H^,. 

Table|5] gives the orbital elements for the close system ob- 
tained by KOREL and ORBITX, respectively. The epoch T 
refers to the time of the maximum positive velocity of compo- 
nent Ba. In the ORBITX analysis, the eccentricity was found to 
be insignificant and was therefore fixed to zero, whereas in the 
KOREL analysis the eccentricity is small, but significant. Fig.|9] 
shows the orbit solution obtained by ORBITX for the narrow- 
line pair (circles) and the velocity of the broad-line component 
(triangles) phased with the ephemeris (HJD = 2453 069.089(3) 
H- 3.5744(3)£') of the close orbit. 

For HD 61199 Ba and Bb, lo wer limi t s on th e masses can be 
derived via the equation given bv lAitkenl(ll935h 



Table 5. Results for the orbit solution of the close binary. 



KOREL 



ORBITX 



P[d] 


3.5720 ± 


0.0015 


3.57436 + 0.00034 


T (2453 000+) 


69.1629 ± 


: 0.0028 


69.0889 ± 0.0025 


e 


0.01003 ± 


0.00036 


(fixed) 




7.59 ± 


0.12 




Ki [kms-'] 


74.564 i 


: 0.038 


74.55 ± 0.24 


Kz [kms-'] 


77.00 i 


:0.18 


77.05 ± 0.27 


Vo [kms-'] 






38.32 ± 0.12 



m\ 2 sin I 



10 



-6.98358 



where m denotes the mass in solar units, K\ and repre- 
sent the semi-amplitudes of the radial velocities in kms-^ The 
period in days is denoted by P and the eccentricity is expressed 
by e. The exponent has not been rounded off intentionally, 
to keep the formula true to original. The lower limits on the 
masses are: for component Ba, 0.68 M©, and for component 
Bb, 0.65 Mq, which, together with similar luminosities deduced 
from line strength, suggests similar values of Tetf and log g for 
both components. 
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Fig. 7. Top panel: Fourier transform of the time-dependent ve- 
locity observable, as explained in the text. The two highest peaks 
correspond to periods of 4.54 and 3.56 days. Note that the radial 
velocity observable varies with twice the orbit frequency. The 
middle and bottom panels show the radial velocity curves of 
the three components of HD61199 phased with 4.54 and 3.56 
days respectively. The filled circles correspond to component 
Ba (which has systematically stronger spectral lines, as seen in 
Figure|6] ); open circles correspond to component Bb; and the 
triangles to the fast rotating component A. 
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Fig. 8. Disentangled spectra of the components Ba (1), Bb (2) , 
A (3, from top to bottom) and the combined spectrum. 
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Fig. 9. Radial velocity curves of the HD61199 system phased 
with the final orbit solution. The filled circles, open circles and 
triangles are the measured radial velocities of the components 
Ba, Bb and A, respectively. The solid and dashed lines show the 
orbital solution of components Ba and Bb in the close binary. 
The phase is based on the ephemeris derived by the ORBITX 
solution. 



Aitken's formula for the masses of binaries allows us to con- 
strain the inclination angle of the close binary orbit. A conserva- 
tive mass range for F type stars ranging from 1 . 1 to 1 .6 yields 
a range for the inclination angle of the orbital plane between 47 
and 57 degrees. 

4.2. Determination of the Fundamental Parameters 
4.2.1. Direct Method 

Further characteristics extracted from the LSD profiles are the 
mean radial velocity for component A of 45.8 ± Ikms"', while 
the 7 velocity of the close binary Ba and Bb was confirmed 
to be 38.2+ Ikms"'. A mean vsin/ was derived for Ba to 
be 15.1 ± l.lkms-i; for Bb, 15.2 ± 1.2kms-'; and for A, 
131 + 1.3 kms '. The relative equivalent width contributions are 
10 + 1 % for Ba, 9 + 1 % for Bb, and 8 1 ± 3% for A. The analysis 
of the Ha region was our main method of finding fundamental 



parameters of the components. The observed spectrum was 
fitted by scaling the individual model spectra using a parameter 
a, defined as follows: 

« = ^i'^L = Rlt = ¥' where + rI„ +R\^1. 

The value of a was adjusted simultaneously with effective 
temperatures. We arrived at a final value of a = 0.7. If the 
system constitutes a physical triple system, the parameter a 
gives the radii ratio RA/Rsa.Bh - VO.7/0.15 = 2.2. 

By comparing fits of synthetic triple-star spectra across a 
wide parameter space, we estimate for component A an effec- 
tive temperature of = (8000 + 200) K and for each compo- 
nent of the close binary Teff - (6400 ± 400) K. The rather large 
uncertainties of the temperatures for the close binary is caused 
by the uncertainty in log g and the fact that its contribution to 
the total light is small compared to the dominant A component. 
Fig. [To] illustrates the quality of the best fit and the accuracy to 
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Fig. 10. Fits to the Ha line of the integrated spectrum of 
HD 61 199. The black lines represent the observations; the solid 
line corresponds to a synthetic spectrum with Tes = 8000 K, 
log g = 4.2, plus two stars with 6800 K (4.4) and 6600 K (4.4). 
Another model (dashed line) with 8000 K (4.0), plus two stars 
with each 6000 K (4.2), also matches the observed spectrum. 
Varying narrow lines are telluric in origin. 



which fundamental parameters of the three components can be 
determined from existing data. The black wiggly lines represent 
four out of 29 observed spectra in the range around Ha and 
the smooth blue line corresponds to a synthetic spectrum with 
Teff = 8000K and log § = 4.2, plus two stars with 6800 K (4.4) 
and 6600 K (4.4). Another model with 8000 K (4.0) plus two 
stars (each with Teff = 6000 K and log g = 4.2) also agrees with 
the observed spectra within the resolution. The radial velocity 
values for these fits were taken from our orbital solution (see 
Table 5). The chosen parameters for the triple system give also 
a good fit to the metal line spectrum as is illustrated in FigJTTI 
Based on evolutionary tracks published by Schaller et al. (1992), 
Fig.fT2]shows that a triple system with log(age) between 8.8 and 
8.9 is compatible with our observations of HD 61199. 

4.2.2. Disentangled Spectra 

The determination of the fundamental parameters was also per- 
formed on the disentangled spectra. A grid of model spectra 
ranging from 6000 to 8600 K with stepwidth of 200 K and log 
g of 3.6, 4.0, 4.4, respectively, was compared to the disentan- 
gled spectra corrected for radial velocity shi fts. The model s pec- 
tra were calculated using the co de Synth V (Tsvmbal"1996') and 
NEMO model atmosphere grid (" Heiter et al. 2002). For finding 
the best fitting spectra two approaches have been used. First, the 
standard deviation of the residuals, hereafter referred to as cr, 
and second, the test were used for evaluating the quality of 
the fit. These approaches lead to different best fitting solutions, 
because the methods differ in weighting the data. The standard 
deviation of the residuals weights each data point equally, while 
the;i'^ test weights the data according to the relative line strength. 



Z = 0.018 




9000 B500 BOOO 7500 7000 6500 6000 



Fig. 12. Evolutionary tracks from ^ Schaller et al.l (11992 ^ for 
masses Mq = 1.25 - 2.0 and various logg values, in the ranges 
we infer for HD61199 based on the spectral fits. The vertical 
dotted lines give the temperatures ranges for the individual com- 
ponents of HD 61 199. The dashed lines correspond to isochrones 
with the given ages at the lower left corner. The isochrone for log 
age of 7.74 corresponds to the lowest line. 



Hence the;i'^ test prefers better fits of line cores. Solar chemical 
composition was assumed, based on the good fit obtained by the 
direct method (see Fig.fTTTi. 

The standard deviation of the residuals and was calcu- 
lated the following way: 



where S seal means scaled synthetic spectrum and D the 
disentangled spectrum and R - S seal - D. The model spectra 
were scaled according to 



where S denotes the model spectrum and / corresponds 
to the luminosity factor This factor is the ratio of the luminosity 
of a given component to the sum of the luminosities of both 
other components. For each model spectrum the best fitting 
luminosity factor / and the optimal v sin / has been determined. 
There are spectral regions where the agreement of synthetic and 
disentangled spectra is rather poor, which may be due to an 
incomplete input line list. Those regions have been omitted for 
finding the final solution. Table|6] gives the best fitting model 
spectra for each component, which are illustrated in Fig. [13] As 
a conclusion, fitting model spectra to the disentangled spectra 
gives about the same result as the direct method discussed 
above. Both fainter components are cooler stars, while the 
brighter star is definitely hotter. The model spectra for the cooler 
stars are only weakly sensitive to changes in logg, while the 
for the hotter component model spectra with lower log g values 
give better cr oix^ values. 

Table[8] shows the results of both analyses, disentangled 
spectra and Hq. fitting of the combined spectrum. The most likely 
parameters are printed bold face. For the fainter components 
large uncertainties in the parameters are evident. 
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Fig. 11. Comparison of our best-fitting 
synthetic spectrum to an observed 
spectrum of HD 61 199 with high S/N 
during a phase of large separation 
of components in the close binary, 
in a wavelength range with many 
metal lines. The model parameters 
for the synthesis are T^ff = 8000 K, 
(log g = 4.2), 6800 K (4.4) and 6600 K 
(4.4), for A, Ba, Bb respectively. 



Table 6. Results for the individual spectra obtained by KOREL. 



Component, method 


Tcff 


V sini 


f 


logg 


A, cr 


8200 


125 


3.7 


3.6 




8000 


131 


8.5 


3.6 


Bb,o- 


6400 


13 


0.04 


4.4 


Bb,^ 


6000 


16 


0.036 


4.4 


Ba, cr 


6400 


13 


0.056 


4.4 




6000 


15 


0.06 


4.4 



5. Asteroseismic modelling of HD61199 A 

Theoretical pulsation frequencies of HD 6 11 99 were com puted 
with updat ed codes for stellar e volution ( |Paczvnskilll969i) and 
oscillation (lDziemblowskill977h. where rotation is tak en into ac- 
count to second order ('Dziemb owskT .&Goode'T992^. Standard 
Mixin g Length Theory (MLT) for convection (Bohm-Vitense 
[T958h was used and the simplest convection - pulsation treat- 
ment (the assumption of frozen convective flux). Since enve- 
lope convection becomes less efficient for stars closer to the hot 
border of the classical instability strip, omlt can be expected 
to be smaller for HD61199 than for the 6 Scuti star FG Vir 
For that star, the MLT parameter o mlt was found to be 0.5 
jPaszvnska-Daszkiewicz et al. I l2005h . 

Direct mode identification based on our photometric and 
spectroscopic data is not possible, so we instead investigate pe- 
riod ratios. The values of radial-pulsation period ratios depend 
mainly on metallicity and rotation (Suarez et al. 2006) and also 
slightly on the evolutionary stage of the star. Table|2] shows 
the predicted values for radial period ratios of a 5 Scuti star 
model with Z = 0.018 and 0.02 and Vrot(ZAMS) = ISOkms ' 
and ISOkms '. 

For convenience, the radial fundamental mode is hereafter 
referred to as F, the first overtone as IH, and so on. The observed 
pulsation frequencies in Table[3]were examined to identify pos- 
sible period ratios for radial modes in HD 61 199. (The two long 
periods given in Table[3] have been excluded because they are 
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Fig. 13. Model spectra compared to the disentangled spectra. 
From top to bottom the comparison is shown for component Ba, 
Bb and A respectively. The full line corresponds to the best fit- 
ting model spectrum using cr, the gray dashed line corresponds 
to the results based on analysis. Details on the fitted model 
spectra are given in Table|6] Notice, that the scaling of the wave- 
length axis is different for A and another wavelength region is 
shown. 



Table 7. Expected period ratios from models. 





Viot =150 kms ' 


Viot = 150 kms ' 


Vi-o, = 180 kms-' 




Z = 0.018 


Z = 0.02 


Z = 0.018 


F/IH 


0.782 


0.782 


0.786 


1H/2H 


0.812 


0.812 


0.813 


2H/3H 


0.842 


0.841 


0.843 


3H/4H 


0.863 


0.862 


0.864 


4H/5H 


0.880 


0.880 


0.882 


5H/6H 


0.892 


0.892 


0.893 
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not due to 6 Scuti pulsation.) Our temperature estimate indicates 
that the star is located close to the hot border of the instability 
strip of radially pulsating 6 Scuti stars (Pamyatnylch 2000). 
Under the assumption that radial modes are observed, we find 
five possibilities: 

(a) /4//5(lH/2H) = 0.819 

(b) /8//„ (2H/3H) = 0.837 

(c) /i2//5(2H/3H) = 0.846 

(d) /2//io (2H/3H) = 0.841 

(e) /„//i (4H/5H) = 0.884 

In f ast rotating stars, nonradial modes usually dominate 
( lBrege i<^2000). and this may also be the case for HD 61199. 
Unfortunately, fitting the observed frequencies to model fre- 
quencies with a least-squares technique does not result in a 
unique solution, because models predict a very dense spectrum 
of frequencies of which only a few (in this case, 11) are detected. 
Hence we concentrated our modelling efforts on explaining the 
observed frequency range of 18.69 - 29.70 1/ '. Figure[T4lshows 
the region in the log g - log Teff diagram where pulsation models 
are consistent with the frequency range observed in HD 61 199. 
The inner polygons are calculated with Z = 0.018 and Z = 0.02 
and represent the area in which all observed frequencies lie in 
the theoretically predicted range of excited modes. Due to uncer- 
tainties in stellar modelling and potentially incomplete opacity 
data the frequency range of unstable modes may be larger than 
predicted by our codes. The outer polygons shows the range of 
models where the observed frequency range is 3 ' larger than 
the theoretically predicted range of unstable modes. 

The models (a) to (e) mentioned above are also included in 
the figure. Models (b), (c) and (d) lie outside the errorbox, model 
(a) is close to the edge of the error box, while (e) is inside. In 
model (e) the period ratio of the fourth overtone and the fifth 
overtone of the radial fundamental would match the theoretical 
ratio given in Table|2] However, the unstable modes range from 
19.6(i"' to 27.3c/"', which is too small by 1 c/"' on the lower 
frequency border and 2.3 c/ ' on the higher frequency border, re- 
spectively. The model parameters are Teff = 8050 K, log g - 3.86, 
M = 2.1 Mq. This model is compatible to the derived parame- 
ters for HD61199 in Sect.|4] The amplitude of the presumably 
radial mode in model (e), /n, is 0.16 mmag and did not change 
significantly, while the amplitude of the presumably consecutive 
radial mode, /i, is about 1.5 mmag and was clearly variable. 
Considering that /i i is close to the detection limit while /i has 
the highest amplitude, it is questionable, if the identification with 
radial modes is correct in this case. 

The error box in Fig. [12] is based on the assumption, that 
HD 61 199 forms a physical triple system, where for component 
A log g =4.0 was adopted. Models hotter than logTeff = 3.92 
cannot explain the observed frequency range, nor can models 
with masses larger than 2.1 Mq. Note, however, that models gen- 
erated for EG Vir tend to be too cool by 200 to 300 K, when 
comparing the best-fitting models to the spectroscopically deter- 
mined temperatures. 

6. Discussion and Conclusions 

HD 61 199 was discovered to be a multiperiodic 6 Scuti pulsator 
from MOST photometry in 2004. Analysis of these data, and 
MOST data collected in 2005 and 2007, reveal 1 1 pulsation fre- 
quencies, with amplitudes as low as 90 ppm in the total flux of 
the system. Subsequent Tautenburg spectroscopy demonstrates 
that HD61199 is a triple system, where the two form a close 
binary with a period of 3.57 days, and the 6 Scuti star (which 
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Fig. 14. Region in the log g - log T^s diagram where models fit 
the observed range of unstable mode frequencies in HD 61 199. 
Evolutionary tracks for M=1.5, 1.75 and 2.0 are given, 
where the full lines coiTespond toZ- 0.018 and the dotted lines 
to Z- 0.02. The vertical dashed line indicates the estimated ef- 
fective temperature of component A, and the error limits of the 
temperature of HD 61 199 A derived in Sec. 14. 21 are indicated 
by the full vertical lines. The filled circles labelled (a) to (e) 
are models which fit possible observed period ratios for radial 
modes (see text). The inner polygon shows the region of models 
where the observed range of unstable frequencies is exactly re- 
produced. The outer polygon shows the region where a deviation 
of 3 d ' is allowed. Both polygons have been calculated with Z 
= 0.018 (marked by the dashed-dotted line, crosses indicate the 
position of calculated models) and Z - 0.02 (marked by dotted 
line and pluses). 



we dub the A component and which contributes 80% of the to- 
tal flux) is either in a much wider orbit or is incidentally in the 
line of sight. The system does not exhibit eclipses, but variations 
with the period of the close binary orbit are evident in the 2007 
MOST photometry. Another period of about 4 days is present in 
all three epochs of photometry. The origin of this period remains 
unknown. The high v sin / of the 6 Scuti star, and its likely ra- 
dius, mean that its rotation period must be about 1 day. The only 
way we know to associate a 4-day period with the 6 Scuti star is 
through an isolated high-overtone g-mode. It could be a hybrid 
pulsator, with both y Doradus g-modes and 6 Scuti p-modes. 
However, the expected range of high order g-mode periods is 
about 0.3 - 3 d, and most y Dor periods are closer to 0.8 d. Thus 
the hybrid explanation is unlikely, and we conclude the 4-day 
period originates in one or both components of the close binary. 
However, the stars in the close binary are almost certainly tidally 
locked, so it is unlikely that the 4-day period represents evidence 
for rotation in that system. If the rotation axes of the two compo- 
nents were significantly misaligned, we would be able to mea- 
sure the difference in v sin / between the two stars in the spectra, 
and we do not see such a difference. 

An additional puzzle in the system is the fact that the orbital 
frequency of the close binary is equal to the difference between 
the two 6 Scuti frequencies with the highest amplitudes. The sep- 
aration of the 6 Scuti star from the close binary is large enough 
that there should be no significant tidal influence. The observed 
frequency coincidence is remarkable, but we have no explana- 
tion for it. 
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Table 8. Summary of fundamental parameters derived with different techniques and proposed solution. 
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Component 




A 






Bb 






Ba 




Method 






V sin i 






V sin i 




log 4' 


V sin i 


KOREL (a) 


8200 


3.6 


125 


6400 


4.4 


13 


6400 


4.4 


13 


KOREL (x^) 


8000 


3.6 


131 


6000 


4.4 


16 


6000 


4.4 


15 


Ua 


8000 


4.0 


131 


6600 


4.4 


15 


6800 


4.4 


15 




8000 


4.2 


131 


6000 


4.2 


15 


6000 


4.2 


15 


Seismology' max 


8100 


4.2 
















min 


7500 


3.85 
















Triple system^ max 




4.3 






4.4 






4.4 




min 




3.9 






4.3 






4.3 




Proposed solution 


8000 


4.0 


130 


6400 


4.4 


15 


6400 


4.4 


15 


estimated uncertainty 


200 


0.5 


3 


400 


0.5 


2 


400 


0.5 


2 



from frequency ranges predicted by theory, see Fig|14l - for chosen temperatures, see Fig.ir2T 



Considering the various estimates for the fundamental astro- 
physical parameters of the three components we conclude that 
(reffAog^) A: 8000/4.0, Ba and Bb: 6400/4.4, are the best com- 
promise within the uncertainties given in Table |8] This conclu- 
sion is also supported by the fact that synthetic spectra of a triple 
system computed for various combinations in the mentioned pa- 
rameter space overlap within the observational uncertainties of 
our spectroscopy. The Hq. line profiles are rather insensitive to 
log g for stars with Tetf < 8500 K, hence the large uncertainty for 
log g given in the table are estimated from the metal lines of the 
disentangled spectra. 

The parallax given bvl Gatewood & Hani (|2006|) is consistent 
with the A component being a 6 Scuti star with an absolute mag- 
nitude of My a: 1.8. For both fainter components, Ba and Bb, 
we derive My ~ 4.1 and 4.2, respectively, based on our model 
of the triple system. S Scuti models enable us to estimate upper 
limits on the effective temperature (about 8000 K) and mass (2.1 
Mq) of HD 61 199 A. Based on the orbit solution, we set lower 
mass limits of about 0.7 Mq for the close binary components, 
Ba and Bb. However, given the expected masses for F-type stars 
(between 1 . 1 and 1 .6 Mq), we can exploit the upper limit to con- 
strain the inclination of the orbital plane to 47° - 57°. 
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